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Abstract 

We study the constraints of the generic two-Higgs-doublet model (2HDM) type-III and the impacts of 
the new Yukawa couplings. For comparisons, we revisit the analysis in the 2HDM type-II. To understand 
the influence of all involving free parameters and to realize their correlations, we employ y-square fitting 
approach by including theoretical and experimental constraints, such as S, T, and U oblique parameters, the 
production of standard model Higgs and its decay to 77, WW*fZZ*, T+r“, etc. The errors of analysis are 
taken at 68%, 95.5%, and 99.7% confidence levels. Due to the new Yukawa couplings being associated with 
cos{P — a) and sin(/3 — a), we find that the allowed regions for sin a and tan (3 in the type-III model can be 
broader when the dictated parameter xf is positive; however, for negative xf^ the limits are stricter than 
those in the type-II model. By using the constrained parameters, we find that the deviation from the SM in 
the h —)■ Z7 can be of 0(10%). Additionally, we also study the top-quark flavor-changing processes induced 
at the tree level in the type-III model and find that when all current experimental data are considered, we 
get Br(t —)• c{h,H)) < 10“^ for mn = 125.36 and mn = 150 GeV and Br{t —)■ cA) slightly exceeds 10“^ for 
TOA = 130 GeV. 
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I. INTRODUCTION 


A scalar boson around 125 GeV was observed in 2012 by ATLAS [I] and CMS | 2 ] at CERN with 
more than 5a significance. The discovery of such particle was based on the analyses of following 
channels: 77 , WW*, ZZ* and with errors of order of 20-30% and bb channel with an error 

of order of 40-50%. The recent updates from ATLAS and CMS with 7 © 8 TeV data [^|7] indicate 
the possible deviations from the standard model (SM) predictions. Although the errors of current 
data are still somewhat large, the new physics signals may become clear in the second run of the 
LHC at 13-14 TeV. 

It is expected that the Higgs couplings to gauge bosons (fermions) at the LHC indeed could 
reach 4-6% (6-13%) accuracy when the collected data are up to the integrated luminosity of 300 
fb-i IHIE]. Furthermore, e^e Linear Collider (LC) would be able to measure the Higgs couplings 
at the percent level m- Therefore, the goals of LHC at run H are (a) to pin down the nature 
of the observed scalar and see if it is the SM Higgs boson or a new scalar boson; (b) to reveal 
the existence of new physics effects, such as the measurement of flavor-changing neutral currents 
(FCNCs) at the top-quark decays, i.e. t —)■ qh. 

Motivated by the observations of the diphoton, WW*, ZZ*, and r+r~ processes at the ATLAS 
and CMS, it is interesting to investigate what sorts of models may naturally be consistent with 
these measurements and what the implications are for other channels, e.g. h —?■ Z 7 and t —)• ch. 
Although many possible extensions of the SM have been discussed mm, it is interesting to 
study the simplest extension from one Higgs doublet to two-Higgs-doublet model (2HDM) [T3l[20] . 
According to the situation of Higgs fields coupling to fermions, the 2HDMs are classified as type-I, 
-H, and -HI models, lepton specific model, and flipped model. The 2HDM type-HI is the case 
where both Higgs doublets couple to all fermions; as a result, FCNCs at the tree level appear. The 
detailed discussions on the 2HDMs are shown elsewhere [El- 

After scalar particle of 125 GeV is discovered, the implications of the observed /i —)• 77 in the 
type-I and H models are studied [ 2 T] and the impacts on L —)• 7 Z are given [ 22 l |23] . As known that 
the tan and angle a are important free parameters in the 2HDMs, where the former is the ratio of 
two vacuum expectation values (VEVs) of Higgses and the latter is the mixing parameter between 
the two CP-even scalars. It is found that the current LHC data put rather severe constraints on 
the free parameters [16]. For instance, the large tan/3 ~ mt/nib scenario in the type-I and -H is 
excluded except if we tune the a parameter to be rather small a < 0.02. Nevertheless, both type-I 
and type-H models can still fit the data in some small regions of tan/3 and a. 
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In this paper, we will explore the influence of new Higgs couplings on the h —)• h —>■ 

gg, 77 , WW, ZZ and h —)• Z 7 decays in the framework of the 2HDM type-III. We will show what 
is the most favored regions of the type-III parameter space when theoretical and experimental 
constraints are considered simultaneously. FCNCs of heavy quark such as t ^ qh have been 
intensively studied both from the experimental and theoretical point of view [2l]. Such processes 
are well established in the SM and are excellent probes for the existence of new physics. In the 
SM and 2HDM type-I and -II, the top-quark FCNCs are generated at one-loop level by charged 
currents and are highly suppressed due to the GIM mechanism. The branching ratio (BR) for 
t —)■ c/i in the SM is estimated to be 3 x 10“^^ [2^. If this decay t ^ ch is observed, it would be 
an indisputable sign of new physics. Since the tree-level FCNCs appear in the type-III model, we 
explore if the Br{t —)• ch) reaches the order of 10“^-10“^ [261127) . the sensitivity which is expected 
by the integrated luminosity of 3000 fb~^. 

The paper is organized as follows. In section II, we introduce the scalar potential and the Yukawa 
interactions in the 2HDM type-III. The theoretical and experimental constraints are described in 
section III. We set up the free parameters and establish the y-sqnare for the best-fit approach 
in section VI. In the same section, we discuss the numerical results when all theoretical and 
experimental constraints are taken into account. The conclusions are given in section V. 


II. MODEL 


In this section we define the scalar potential and the Yukawa sector in the 2HDM type-III. The 
scalar potential in SU{2)l ® U{1)y gauge symmetry and CP invariance is given by [32] 


V($l,$2) 


= + h.c) 

-h ^A2($2^2)^ + A3 ($|^>i)(^>^$2) + A4($1^>2)(^1^2) 










-|- A7^>2^2 


<I>|<I >2 + h.c. 


( 1 ) 


where the doublets <ki ^2 have weak hypercharge Y = 1, the corresponding VEVs are vi and ^ 2 , and 
Aj and m \2 are real parameters. After electroweak symmetry breaking, three of the eight degrees 
of freedom in the two Higgs doublets are the Goldstone bosons (G^, and the remaining five 
degrees of freedom become the physical Higgs bosons: 2 CP-even h, H, one CP-odd A, and a pair of 
charged Higgs H^. After using the minimized conditions and the W mass, the potential in Eq. Q 
has nine parameters which will be taken as: (Ai)j=i^,,, 7 , 771 ^ 2 , and tan/? = V 2 /v\. Equivalently, we 
can use the masses as the independent parameters; therefore, the set of free parameters can be 
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chosen to be 


{ruh , ruH , mA , m^± , tan /3 , a , 771 ^ 2 } > 


( 2 ) 


where we only list seven of the nine parameters, the angle 13 diagonalizes the squared mass matrices 
of CP-odd and charged scalars and the angle a diagonalizes the CP-even squared-mass matrix. In 
order to avoid generating spontaneous CP violation, we further require 

Aeuf 


mi2 - 


hvl 


> C>^^ViV2 


( 3 ) 


with C = 1(0) for A 5 > (<)0 [32]. It has been known that by assuming neutral flavor conservation 
at the tree-level [IHIj we have four types of Higgs couplings to the fermions. In the 2HDM type-I, 
the quarks and leptons couple only to one of the two Higgs doublets and the case is the same as the 
SM. In the 2HDM type-H, the charged leptons and down type quarks couple to one Higgs doublet 
and the up type quarks couple to the other. The lepton-specific model is similar to type-I, but the 
leptons couple to the other Higgs doublet. In the flipped model, which is similar to type-H, the 
leptons and up type quarks couple to the same double. 

If the tree-level FCNCs are allowed, both doublets can couple to leptons and quarks and the 
associated model is called 2HDM type-HI [T5l|29l[30]. Thus, the Yukawa interactions for quarks 
are written as 


= QLY^dR(j)k + QiY^UR^k + h.c. (4) 

where the flavor indices are suppressed, = {ur, dL) is the left-handed quark doublet, and 
denote the 3x3 Yukawa matrices, = ia 2 <p%, and k is the doublet number. Similar formulas could 
be applied to the lepton sector. Since the mass matrices of quarks are combined by Y^{Y^) and 
y 2 (y 2 ) fjown (up) type quarks and generally cannot be diagonalized simultaneously, 

as a result, the tree-level FCNCs appear and the effects lead to the oscillations of K — K, Bq — Bq 
and D — D at the tree-level. To get naturally small FCNCs, one can use the ansatz formulated 
by yA, Yjj oc ^rriimj/v |29l[30|. After spontaneous symmetry breaking, the scalar couplings to 
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fermions can be expressed as m 


--2HDM-III 

l-y 


= ULi 


COS a mu, . cos(^ - «) i u 


sin /3 V 


-\/2sin /3 




cos (3 V 


y/2 cos 13 


+ ULi 
+ dii 
- iULi 


sin a mu, ^ , sin(/3 -a) 

cosamd sin(/3-a) ^ ^ 

1 m,,. 




^ _ 

tan/3 u \/2sin/3 


, -1 ( . o^di. , ^ij 

+ idLi -tan/3 —-dij + ^ 

\ V ^J2 cos (3 


URjA 

dRjA + h.c, 


(5) 


where v = ■sjv\^ v\^ X^- = y/n\imq~/vXij {q = u^d ) and xfj are the free parameters. By above 
formulation, if the FCNC effects are ignored, the results are returned to the case of the 2HDM 
type-II, given by 


^2HDM-II 

‘-'Y 


ULi 


f cosamuj ^ \ 

Vsin,0 V *■^7 


URj h “I" dRi 


sin a mdi 
cos /3 V 



dRjh + h.c. 


( 6 ) 


For the couplings of other scalars to fermions, the can be found elsewhere |31) . It can be 
seen clearly that if x^j'^ are of 0(10“^), the new effects are dominated by heavy fermions and 
comparable with those in the type-II model. The couplings of h and H to gauge bosons V = W, Z 
are proportional to sin(/3 — a) and cos(/3 — a), respectively. Therefore, the SM-like Higgs boson 
h is recovered when cos(/3 — a) ~ 0. The decoupling limit can be achieved if cos(/3 — a) ~ 0 
and mh 'C mH,mA,mH± are satisfied [32]. From Eqs. (§ and Q, one can also find that in the 
decoupling limit, the h couplings to quarks are returned to the SM case. 

In this analysis, since we take a in the range —7r/2 < a < tiI2, sin a will have both positive and 
negative sign. In the 2HDM type-II, if sin a < 0 then the Higgs couplings to up- and down-type 
quarks will have the same sign as those in the SM. It is worthy to mention that sin a in minimal 
supersymmetric SM (MSSM) is negative unless some extremely large radiative corrections flip its 
sign [32]. If sin a is positive, then the Higgs coupling to down quarks will have a different sign with 
respect to the SM case. This is called by the wrong sign Yukawa coupling in the literature [321133] . 
Later we will explain if the type-HI model would favor such wrong sign scenario or not. 
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III. THEORETICAL AND EXPERIMENTAL CONSTRAINTS 


The free parameters in the scalar potential defined in Eq. 0 could be constrained by theoretical 
requirements and the experimental measurements, where the former mainly includes tree level 
unitarity and vacuum stability when the electroweak symmetry is broken spontaneously. Since the 
unitarity constraint involves a variety of scattering processes, we adopt the results [MIES]. We also 
force the potential to be perturbative by requiring that all quartic couplings of the scalar potential 
obey |Aj| < Svr for all i. For the vacuum stability conditions which ensure that the potential is 
bounded from below, we require that the parameters satisfy the conditions as IMIEZ] 

Ai > 0 , A2 > 0 , As + \/ A1A2 > 0, \/A1A2 + As + A4 — IA5I > 0, 

2\\q + A 7 I < -(Ai + A 2 ) + As + A 4 + As . (7) 

In the following we state the constraints from experimental data. The new neutral and charged 
scalar bosons in 2HDM will affect the self-energy of W and Z bosons through the loop effects. 
Therefore, the involved parameters could be constrained by the precision measurements of the 
oblique parameters, denoted by S, T, and U |38|. Taking mh = 125 GeV, mt = 173.3 GeV and 
assuming that U = 0, the tolerated ranges for S and T are found by |M| 


A5' = 0.06 ±0.09, AT = 0.10 ±0.07, (8) 

where the correlation factor is p = ±0.91, A5 = 5 ' 2 HDM _ ^SM _ ji 2 HDM _ j^SM^ their 

explicit expressions can be found [32l • We note that in the limit mfj± = m^o or mH± = rrijjo , AT 
vanishes jiniill]. 

The second set of constraints comes from B physics observables. It has been shown recently in 
Ref. |l2| that Br{B —>■ gives a lower limit on mjj± > 480 GeV in the type-II at 95%GL. By 

the precision measurements of Z ^ bb and Bq — Bq mixing, the values of tan/3 < 0.5 have been 
excluded |l3|. In this work we allow tan^S > 0.5. Except some specific scenarios, tan/3 can not be 
too large due to the requirement of perturbation theory. 

By the observation of scalar boson at mh ~ 125 GeV, the searches for Higgs boson at ATLAS and 
CMS can give strong bounds on the free parameters. The signal events in the Higgs measurements 
are represented by the signal strength, which is defined by the ratio of Higgs signal to the SM 
prediction and given by 


/ _ (Ti{h) • Br{h /) 

- af^{h)-BrSMi^h^f) 


— ■ Pf 1 


(9) 
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where ai{h) denotes the Higgs production cross section by channel i and Br{h —)■ /) is the BR 
for the Higgs decay h ^ f. Since several Higgs boson production channels are available at the 
LHC, we are interested in the gluon fusion production {ggF), tth, vector boson fusion (VBF) and 
Higgs-strahlung Vh with V = WjZ; and they are grouped to be and ^^ypp^vh- order 

to consider the constraints from the current LHC data, the scaling factors which show the Higgs 
coupling deviations from the SM are defined as 


„2HDM 
_ yhVV 

Ky = Kw = KZ = 4 ; ■ 

9hVV 


Kf = 


„,2HDM 

_ yhff 

„SM • 

yhff 


( 10 ) 


where ghw and yhff are the Higgs couplings to gauge bosons and fermions, respectively, and / 
stands for top, bottom quarks, and tau lepton. The scaling factors for loop-induced channels are 
defined by 


III 

(M 

r(h — 77 ) 2 hdm 


r(h —)• g g) 2 HDM 

r(/i-^77)sM 

^{h gg)sM 

2 

T{h ^7)2HDM 

2 

r(/i)2HDM 

— 




T{h -)• 2 ' 7 )sm 

r(^)sM 


( 11 ) 


where r(/i — )• XY) is the partial decay rate for h —?■ XY. In this study, the partial decay width of 
the Higgs is taken from [H], where QCD corrections have been taken into account. In the decay 
modes h —)• 77 and h —?• Z 7 , we have included the contributions of charged Higgs and new Yukawa 
couplings. Accordingly, the ratio of cross section to the SM prediction for the production channels 
ggF -|- tih and XBF+Vh can be expressed as 


^ ggF+tth 
FyBF+Vh 


i^asMiggP) + K^asMjtth) 

(^sM{ggF) + asM{tth) 

fI^(Tsm{VBF) + Kzh^SM{Zh) + K^asM{Zh) + nl^asMiWh) 
<ySM{VBF) + asM{Zh) + asM{Zh) + crsM{Wh) 


( 12 ) 

(13) 


where asM{Zh) is from the coupling of ZZh and occurs at the tree level and asM{Zh) = asuigg 
Zh) represents the effects of top-quark loop. With mh = 125.36 GeV, the scalar factor Uzh can be 
written as [3] 


Kzh = 2.27K\ + f).2,7Kl-IMKznt. 


(14) 


In the numerical estimations, we use ruh = 125.36 GeV which is from LHG Higgs Cross Section 
Working Group [ 8 j at y/s = 8 TeV. The experimental values of signal strengths are shown in 
Table. where the results of ATLAS [5] and CMS m are combined and denoted by V-ggp^^^h 

y-VEF+Vh- 
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TABLE I: Measured signal strengths /iggF+tth and /tvBF+vh that combine the best fit of ATLAS and CMS 
and correlation coefficient p for the Higgs decay mode [SIIIS]. 


/ 

l^^F+tth 

t^VEF+Vh 

i l<^ggF+tth 

± ICTVBF+Vh 

P 

77 

1.32 

0.8 

0.38 

0.7 

-0.30 


1.70 

0.3 

0.4 

1.20 

-0.59 

WW* 

0.98 

1.28 

0.28 

0.55 

-0.20 

TT 

2 

1.24 

1.50 

0.59 

-0.42 

hi 

1.11 

0.92 

0.65 

0.38 

0 


IV. PARAMETER SETTING, GLOBAL FITTING, AND NUMERICAL RESULTS 


A. Parameters and global fitting 


After introducing the scaling factors for displaying the new physics in various channels, in the 
following we show the explicit relations with the free parameters in the type-III model. By the 


definitions in Eq. (10), the scaling factors for Ky and k/ in the type-III are given by 


Ky = sin(/3 — a ), 


cos a cos(/3 — a) 

KJJ = Kt = Kc= yryns “ XF 


Av Y) — — Av-T” — 


sin /3 
sin a 


+ XF 


y/2 sin 13 ’ 
cos(/3 — a) 


(15) 


cos/3 ' V^cos/3 

f 

Although FCNC processes give strict constraints on flavor changing couplings xij with i ^ j, 
however, the constraints are applied to the flavor changing processes in K, D and B meson systems. 
Since the couplings of scalars to the light quarks have been suppressed by niqjv, the direct limit on 
flavor-conserved coupling xii is mild. Additionally, since the signals for top-quark flavor changing 
processes haven’t been observed yet, the direct constraint on = ^fntmqjvx^i is from the 


experimental bound of t —?■ hqi. Hence, for simplifying the numerical analysis, in Eq. (15) we have 
set X 22 — X 33 = X 33 = X 33 = Xf- Since A 33 = mt/vxF, it is conservative to adopt the vale of 
Xf to be 0(1). In the 2HDM, the charged Higgs will also contribute to h —)> 77 decay and the 
associated scalar triplet coupling hH'^H~ is read by 




1 


2m^ 


cos(a + (3) 
sin 2,0 


{2m\ — 2X^v^) — sin (/0 — a){m\ — 2 m^±) 


^6 


Ay 


sin^ (3 cos^ (3 ^ 


-b cos(/3 — a) 


(16) 
























The scaling factors for loop induced processes h —)• ( 77 , Z 7 , gg) can be expressed by 








1.268kvf ~ 0.279^4 + 0.0042k;, + 0.0034kc + O.OOSGkt — 0.0014A;ij7±j77 

.2 

1.058kvi/ — 0.059k; + O.OOOSGk^ + 0.00014kc — 0.00054A/jJ7±j7=f 

.2 

1.078Ki — 0.065Kb ~ 0.013 kc 


(17) 


where we have used nih = 125.36 GeV and taken m^± = 480 GeV. It is clear that the charged 
Higgs contribution to /i —)• 77 and h — )• Z7 is small. In order to study the influence of new free 
parameters and to understand their correlations, we perform the x-square fitting by using the LHG 
data for Higgs searches mmmm - For a given channel / = 77, WW*, ZZ *, tt , we define the Xf 
as 

1 , f 1 , f 2/3 








■(/^I-Ai)(i^ 2 -^ 2 )- (18) 


aj{l-p‘^) (Ti(T 2(1 - p^) 

where ^{ 2 ^ ^ 1,2 and p are the measured Higgs signal strengths, their one-sigma errors, and their 
correlation, respectively and their values could refer to Table the indices 1 and 2 in turn stand 
for ggF -|- tth and VBF -|- Vh, and 2 are the results in the 2HDM. The global y-square is defined 
by 


2 2 , 2 

X — 2_^Xf + xsT 

f 


(19) 


where the Xst related to the x^ for S and T parameters, the definition can be obtained from 


Eq.(18) by replacing p( —)■ 52 HDM j^ 2 HDM ^ corresponding values can be found 

from Eq. ([^. We do not include bb channel in our analysis because the errors of data are still large. 

In order to display the allowed regions for the parameters, we show the best fit at 68 %, 95.5%, 
and 99.7% confidence levels (GLs), that is, the corresponding errors of x^ are Ax^ < 2.3, 5.99, and 
11.8, respectively. For comparing with LHC data, we require the calculated results in agreement 
with those shown in ATLAS Fig. 3 of Ref. [3] and in CMS Fig. 5 of Ref. [7]. 


V. NUMERICAL RESULTS 

In the following we present the limits of current LHC data based on the three kinds of CL 
introduced in last section. In our numerical calculations, we set the mass of SM Higgs to be 
rrih = 125.36 GeV, and scan the involved parameters in the chosen regions as 

480 GeV < < 1 TeV, 126 GeV < mu < 1 TeV, 100 GeV < niA < 1 TeV , 

-l<sina<l, 0.5 < tan/3 < 50, -(1000GeV)^ < m ?2 < (1000GeV)^ (20) 


9 












The main difference in the scalar potential between type-II and type-III is that the Agj terms 
appear in the type-III model. With the introduction of Agy terms in the potential, not only the 
mass relations of scalar bosons are modified but also the scalar triple and quartic coupling receives 
contributions from Ag and Ay. Since the masses of scalar bosons are regarded as free parameters, 
the relevant Agy effects in this study enters game through the triple coupling that 

contributes to the /i —)• 77 decay, as shown in Eq. (16) and the first line of Eq. Since the 

contribution of the charged Higgs loop to the /i —)• 77 decay is small, expectably the influence of 
Agy on the parameter constraint is not significant. To demonstrate that the contributions of Agy 
are not very important, we present the allowed ranges for tan/3 and sin a by scanning Ag,? in the 
region of [— 1 , 1 ] in Eig.[^ where the theoretical and experimental constraints mentioned earlier are 
included and the plots from left to right in turn stand for = 11.8, 5.99, and 2.3, respectively. 


Additionally, to understand the influence of xf defined in Eq. (15), we also scan xf = [—1,1] in 
the plots. By comparing the results with the case of Agy <C 1 and xf = 1 which is displayed in 
the third plot of Eig. (|^, it can be seen that only small region in the positive sin a is modified and 
the modifications happen only in the large errors of the plot with Ay^ = 2.3 has almost no 
change. Therefore, to simplify the numerical analysis and to reduce the scanned parameters, it is 
reasonable in this study to assume Xqj <C 1. Since the influence of |yir| < 1 should be smaller, to 
get the typical contributions from ECNC effects, we illustrate our studies by setting xf = ±1 in 
the whole numerical analysis. 



sina sina sina 


FIG. 1: The allowed regions in (sin a, tan/3) constrained by theoretical and current experimental inputs, 
where we have used m/j = 125.36 GeV in the type-III with —I < xf < 1 and —I < Agy < 1. The errors for 
y-square fit are 99.7% GL (left panel), 95.5% GL (middel panel) and 68% CL(right panel). 


With Agy <C 1, we present the allowed regions for sin a and tan/3 in Eig. where the left, 
middle and right panels stand for the 2HDM type-II, type-III with xf = ~1 and type-III with 
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Xf = +1, respectively, and in each plot we show the constraints at 68% CL (green), 95.5% CL (red) 
and 99.7% CL (black). Our results in type-II are consistent with those obtained by the authors 
in Refs. HSIEI] when the same conditions are chosen. By the plots, we see that in type-III with 
Xf = — 1) ) due to the sign of coupling being the same as type-II, the allowed values for sin a and 
tan/3 are further shrunk; especially sin a is limited to be less than 0.1. On the contrary, type-III 
with Xf = +1) the allowed values of sin a and tan/3 are broad. 

As discussed before, the decoupling limit occurs at a —>■ /3 —f/2 , i.e. sin a = — cos/3 < 0. Since 
we regard the masses of new scalars as free parameters and scan them in the regions shown in 
Eq. (20), therefore, the three plots in Fig. cover lower and heavier mass of charged Higgs. We 


further check that sin a > 0 could be excluded at 95.5(99.7)% CL when > 585(690) GeV. 
The main differences between type-II and type-III are the Yukawa couplings as shown in Eq. ([^ . 

50 r 



sina 




- 0.5 0.0 0.5 

sinor 


FIG. 2: The allowed regions in (sin a, tan/3) constrained by theoretical and current experimental inputs, 
where we have used ruh = 125.36 GeV, the left, middle and right panels stand for the 2HDM type-II, type-III 
with Xf = —1 and type-III with xf = +1, respectively. The errors for y-square fit are 99.7% CL (black), 
95.5% CL (red) and 68% CL (green). 


In order to see the influence of the new effects of type-III, we plot the allowed Kg as a function 
of sin a and tan /3 in Fig. where the three plots from left to right correspond to type-II, type-III 
with Xf = —1 and type-III with xf = +1, the solid, dashed and dotted lines in each plot stand for 
the decoupling limit (DL) of SM, 15% deviation from DL and 20% deviation from DL, respectively. 
For comparisons, we also put the results of 99.7% in Fig. [^in each plot. By the analysis, we see 
that the deviations of Kg from DL in xf = +1 are clear and significant while the influence of 
Xf = — 1 is small. It is pointed out that a wrong sign Yukawa coupling to down type quarks could 
happen in type-II 2HDM [321133) . For understanding the sign flip, we rewrite the kd defined in 
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FIG. 3: Kg as a function of sin a and tan/3 in type-II (left) and type-III with xf = (—lj+1) (middle, 
right), where the solid, dashed and dotted lines in each plot stand for the decoupling limit (DL) of SM, 15% 
deviation from DL and 20% deviation from DL, respectively. The dotted points are the allowed values of 
parameters resulted from Fig. 


Eq. (15) to be 


sin a I Xf sin I3\ xf cos a 

fd = -^ 1 1 -) + 


( 21 ) 


COS;0 V* \/2 y \/2 ' 

In the type-II case, we know that in the decoupling limit kd = 1, but kd < 0 if sin a > 0. 
According to the results in the left panel of Fig. sin a > 0 is allowed when the errors of best 


fit are taken by 2 f 7 or 3 ct. The situation in type-III is more complicated. From Eq. (21), we see 


that the factor in the brackets is always positive, therefore, the sign of the first term should be 
the same as that in type-II case. However, due to a G [—'k/2, f/2], the sign of the second term in 


Eq. (21) depends on the sign of xf- For xf = — I, even sin a kd could be negative when the 


first term is smaller than the second term. For XF = +1; if sin a > 0 and the hrst term is over 
the second term, kd < 0 is still allowed. In order to understand the available values of kd when 
the constraints are considered, we present the correlation of kjj and in Fig.[^ where the panels 
from left to right stand for type-II, type-III with xf = and type-III with xf = +1- In each 
plot, the results obtained by x-square htting are applied. The similar correlation of Ky and kd 
is presented in Fig. By these results, we find that comparing with type-II model, the negative 
Kd gets more strict limit in type-III, although a wider parameter space for sin a > 0 is allowed in 
type-III with xf = +1- 

Besides the scaling factors of tree level Higgs decays, kd,u and Ky, it is also interesting to 
understand the allowed values for loop induced processes in 2HDM, e.g. h —)> 77 , gg, and Z 7 , etc. 
It is known that the differences in the associated couplings between /i —)• 77 and gg are the colorless 
W-, T-, and H=*=-loop. By Eq. ( [l7| , we see that the contributions of r and are small, therefore. 
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Kd Kd Kd 


FIG. 4: Correlation of kd and ku, where the left, middle and right panels represent the allowed values 
in type-II, type-III with xf = —1 and type-III with xf = +1, respectively and the results of Fig. are 
applied. 



Kd Kd Kd 

FIG. 5: The legend is the same as that in Fig. but for the correlation of kv and kd- 


the main difference is from the 14^-loop in which the Ky involves. By using the y-square fitting 
approach and with the inputs of the experimental data and theoretical constraints, the allowed 
regions of and Kg in type-II and type-III are displayed in Fig. where the panels from left to 
right are type-II, type-III with xf = — 1 and -|-1; the green, red and black colors in each plot stand 
for the 68%, 95.5% and 99.7% CL, respectively. We find that except a slightly lower k^ is allowed 
in type-II, the hrst two plots have similar results. The situation can be understood from Figs. 
and where the kjj in both models is similar while Ky in type-II could be smaller in the region of 
negative kd', that is, a smaller Ky will lead a smaller k^. In xf = +1 case, the allowed values of 
K^ and Kg are localized in a wider region. 

It is known that except the different gauge couplings, the loop diagrams for h —)■ and 
h ^ XI ^re exact the same. One can understand the loop effects by the numerical form of 
Eq. (|17[). Therefore, we expect the correlation between kz-^ and k^ should behave like a linear 
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Ky Ky Ky 


FIG. 6 : Correlation of and Kg, where the left, middle and right panels represent the allowed values in 
type-II, type-III with xf = —1 and type-III with xf = +1, respectively and the results in Fig. [^obtained 
by x-square fitting are applied. 


relation. We present the correlation between and in Fig.[^ where the legend is the same as 
that for Fig. 1^ From the plots, we see that in most region is less than the SM prediction. The 
type-III with xf = gets more strict constraint and the change is within 10%. For xf = +I 5 the 
deviation of KZ'y from unity could be over 10%. From run I data, the LHC has an upper bounds 
on h —)• Z 7 , at run II this decay mode will be probed. We give the predictions at 13 TeV LHC for 
the signal strength and in Fig. Hence, with the theoretical and experimental 

constraints, is bounded and could be 0(10%) away from SM at 68 %CL. 



Ky Ky Ky 


FIG. 7: The allowed regions in (k^, kzj) plan after imposing theoretical and experimental constrains. Color 
coding the same as Fig. 
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FIG. 8: Correlation between and after imposing theoretical and ex¬ 

perimental constrains. Left, middle and right panels represent the allowed values in type-II, type-III with 
Xf = ~1 and type-III with xf = +1; respectively and the results in Fig. [^obtained by y-square fitting are 
applied. 

VI. t ch DECAY 


In this section, we study the flavor changing t ^ ch process in type-III model. Experimentally, 
there have been intensive activities to explore the top FCNCs. CDF, DO and LEPII collaborations 
have reported some bounds on top FCNCs. At the LHC with rather large top cross section, ATLAS 
and CMS search for top FCNCs and put a limit on the branching fraction which is Br{t —>■ ch) < 
0.82 % for ATLAS |26] and Br{t —> ch) < 0.56 % for CMS [27]. Note that CMS limit is slightly 
better than ATLAS limit. CMS search for t —)• c/i in different channels: h —>■ 77, WW*, ZZ * , and 
T~^T~ while ATLAS used only diphoton channel. With the high luminosity option of the LHC, the 
above limit will be improved to reach about Br{t —?• ch) < 1.5 x lO”'^ [26] for ATLAS detector. 

From the Yukawa couplings in Eq, (§, the partial width for t ^ ch decay is given by 

X ^l-{xh- Xc)^^yl - (Xh + Xc)^ (22) 

where Xc = rric/mt, Xh = mh/mt and is a free parameter and dictates the FCNC effect. It 
is clear from the above expression that the partial width of t —)• c/i is proportional to cos(/3 — a). 
As seen in the previous section, in the case where h is SM-like, cos(/3 — a) is constrained by LHC 
data to be rather small and the t ^ ch branching fraction is limited. As we will see later in 2HDM 
type-H with flavor conservation the rate for t ^ ch is much smaller than type-III |47j . Since we 
assume that the charged Higgs is heavier than 400 GeV, the total decay width of the top contains 
only t ^ ch and t —)■ bW. With nih = 125.36 GeV, mt = 173.3 GeV and rric = 1.42 GeV, the total 
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width can be written as 


r. = n 


SM ^ amn ( 


sin /3 


GeV 


(23) 


where is the partial decay rate for t —)• Wb and is given by 


pSM ^ Gpru^ ^ 2”^^ 


87r\/2 


mt 


m. 


J(l_22fM(2^-()) = l,43GeV 


Stt 


in which the QCD corrections have been included. By the above numerical expressions together 
with the current limit from ATLAS and CMS, the limit on the tch FCNC coupling is found by 
'cos(/3 - a)X^^' 


sin/3 

cos(/3 - a)X |3 
sin/3 


< 2.2 

< 0.36 


for 

for 


Br{t —^ ch) < 8.2 X 10 ^ , 
Br{t —)■ ch) < 5.6 x 10“^ 


(24) 


in agreement with 


We perform a systematic scan over 2HDM parameters, as depicted in Eq. (20), taking into 
account LHC and theoretical constraints. Although X 23 is a free parameter, in order to suppress 
the FCNC effects naturally, as stated earlier we adopt X 23 = vx 23 - Since the current 


experimental measurements only give a upper limit on t —?■ he, basically X 23 is limited by Eq. (24) 
and could be as large as 0(1 — 10^), depending on the allowed value of cos(/3 — a). In order to use 
the constrained results which are obtained from the Higgs measurements and the self-consistent 
parametrisation A 33 = mtjvxF which was used before, we assume X 23 = Xf = il. In our 
numerical analysis, the results under the assumption should be conservative. In Fig. [^left) we 
illustrate the branching fraction of t —)• c/i in 2HDM-HI as a function of cos(/3 — a). The LHC 
constraints within la restrict cos(/3 —a) to be in the range [—0.27,0.27]. The branching fraction for 
t ^ ch is slightly above 10“^ . The actual CMS and ATLAS constraint on Br(t —)> ch) < 5.6 x 10“^ 
does not restrict cos(/3 — a). The expected limit from ATLAS detector with high luminosity 3000 
fb“^ is depicted as the dashed horizontal line. As one can see, the expected ATLAS limit is somehow 
similar to LHC constraints within Icr. In the right panel, we show the allowed parameters space 
in (sina, tan/3) plan where we apply ATLAS expected limit Br{t —)• ch) < 1.5 x 10“^. This plot 
should be compared to the right panel of FigJ^ It is then clear that this additional constraint only 
act on the 3a allowed region from LHC data. 

In Fig. 10 and 11, we show the fitted branching fractions for t ^ ch (left), t —?• cH at 


rnn = 150 GeV (middle) and t —>■ cA at niA = 130 GeV (right) as a function of ku, where Fig. 10 


is for Xf = +1 while Fig. is xf = —1- In the case of xf = +1 the fitted value for kjj at the 3a 
level is in the range [0.6,1.18] and the branching fraction for t —>• ch,cH are less than 10“^ while 


16 

















FIG. 9: Left: Branching ratio of Br(t —>■ ch) as a function of cos(/3 —a), the two horizontal lines correspond to 
LHC actual limit (upper line) and expected limit from ATLAS with 3000 fb“^ luminosity (dashed line). Right 
panel: allowed parameters space in type III with the ATLAS expected limit on Br{t —)■ ch) < 1.5 x 10“^. 

Br{t —)> cA) slightly exceed the 10“^ level. Similarly, for xf = the htted value for ku at the 
3(7 level is in the range [0.85,1.25] and the branching fraction for t —)• ch, cH, cA are the same size 
as in the previous case. 



% % % 


FIG. 10: Branching ratios of Br(t —>■ cft.)(left), Br(t —>■ ciL)(middle) and Br(t —c/i)(right) as a function of 
Kjj in type III with xf = +1- For ruh = 125.36 GeV, mn = 150 GeV and = 130 GeV. 


VII. CONCLUSIONS 

For studying the constraints of 8 TeV LHC experimental data, we perform y-square analysis 
to hnd the most favorable regions for the free parameters in the two-Higgs-doublet models. For 
comparisons, we focus on the type-II and type-III models, in which the latter model not only affects 
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FIG. 11: The legend is the same as Fig.[^but for xf = ~1- 

the flavor conserving Yukawa couplings, but also generates the scalar-mediated flavor-changing 
neutral currents at the tree level. 

Although the difference between type-II and type-III is the Yukawa sector, however, since the 
new Yukawa couplings in type-III are associated with cos(/3 — a) and sin(/3 — a), the modified 
couplings tth and btH^ will change the constraint of free parameters. 

In order to present the influence of modified Yukawa couplings, we show the allowed values of 
sin a and tan/3 in Fig. where the LHC updated data for pp —)• h —;> / with / = 77, WW* jZZ* 
and T~^T~ are applied and other bounds are also included. By the results, we see that sin a and 
tan /3 in type-III gets even stronger constraint if the dictated parameter xf = — 1 is adopted; on the 
contrary, if we take xf = +1; the allowed values for sin a and tan/3 are wider. It has been pointed 
out that there exist the wrong sign Yukawa couplings to down type quarks in the type-II model, 
i.e. sin a > 0 or «£> < 0. By the study, we find that except the allowed regions of parameters are 
shrunk slightly, the situation in xf = is similar to the type-II case. In xf = +1; although the 
< 0 is not excluded completely yet, but the case has a strict limit by current data. We show 
the analyses in Figs.|^and[^ In these figures, one can also see the correlations with modified Higgs 
coupling to top-quark kjj and to gauge boson Ky- 

When the parameters are bounded by the observed channels, we show the influence on the 
unobserved channel h —)■ Z 7 by using the scaling factor KZ'y, which is defined by the ratio of decay 
rate to the SM prediction. We find that the change of KZ'y in type-III with xf = — 1 is less than 
10%; however, with xf = +1) the value of KZ'y could be lower from SM prediction by over 10%. 
We also show our predictions for signal strengths H'y'y and fj^yz and their correlation at 13 TeV. 
The main difference between type-II and -HI model is: the flavor changing neutral currents in 
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the former are only induced by loops, while in the latter they could occur at the tree level. We 
study the scalar-mediated t —?• c{h, H, A) decays in type-III model and find that when all current 
experimental constraints are considered, Br{t —?• c{h,H)) < 10“^ for nifi = 125.36 and rriH = 150 
GeV and Br{t —>■ cA) slightly exceeds 10“^ for niA = 130 GeV. The detailed numerical analyses 
are shown in Figs. [9, 10 and 11 
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